Einstein-Podolsky-Rosen correlations via dissociation of a molecular Bose-Einstein 

condensate 
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Recent experimental measurements of atomic intensity correlations through atom shot noise sug- 
gest that atomic quadrature phase correlations may soon be measured with a similar precision. 
We propose a test of local realism with mesoscopic numbers of massive particles based on such 
measurements. Using dissociation of a Bose-Einstein condensate of diatomic molecules into bosonic 
atoms, we demonstrate that strongly entangled atomic beams may be produced which possess 
Einstein-Podolsky-Rosen (EPR) correlations in field quadratures, in direct analogy to the position 
and momentum correlations originally considered by EPR. 
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The recent demonstrations of atomic correlation mea- 
surements at the shot noise level 0, are a significant 
step towards true quantum atom optics. Quantum op- 
tics, which began with photon correlation measurements, 
has allowed for many fundamental tests of quantum me- 
chanics. Importantly, the availability of lasers allowed 
the development of techniques to perform quadrature 
phase measurements. In quantum atom optics, Bose- 
Einstein condensates (BEC) play the role of the laser. 
However, homodyne noise correlation measurements of 
atomic field quadratures have not been available. 

In this Letter, we suggest one route to achieve this, 
and outline a scheme which would allow for fundamental 
tests of quantum mechanics with massive particles. We 
base our proposal on the Einstein-Podolsky-Rosen (EPR) 
paradox |3j. The EPR paper introduced two particles 
with perfect correlations (entanglement) in momenta and 
positions, these persisting with spatial separation. De- 
pending on which property of one particle we choose to 
measure, we can predict with certainty the same observ- 
able of the other particle. EPR concluded that local real- 
ism was inconsistent with the completeness of the quan- 
tum mechanical description of nature. As suggested by 
Reid in 1989 0], products of variances of inferred optical 
quadratures can demonstrate the paradox by seeming to 
violate the Heisenberg uncertainty relation, although this 
is impossible for directly observed quadratures. This is 
applicable to realistic correlations, and was demonstrated 
experimentally by Ou et al. Q in 1992. 

We show here that dissociation 0, 0, 0, of a molec- 
ular BEC can also exhibit EPR correlations in atomic 
quadratures. Such tests of quantum mechanics (see also 
Ref. [lf]j| ) are a step toward understanding the prop- 
erties of mesoscopic superpositions of massive particles, 
since they introduce couplings to gravitational fields, not 
previously tested in quantum measurement experiments. 

There has been much experimental progress [ill ll2T | 



and intense theoretical interest [T^ . fl^. fl5 l fTH . fl7l TTM 

in the production of molecular dimers 
from Bose condensed atoms. We will assume that this 
can be created from bosonic constituents and propose a 
realization of the EPR paradox via dissociation, which 
can automatically yield two counter-propagating beams 
through momentum conservation. Starting from a three- 
dimensional (3D) molecular condensate, Diirr et al. 
have used Feshbach resonance techniques to create a 
quasi- mono-energetic expanding wave of bosonic atoms, 
in close analogy to successful fermionic correlation ex- 
periments 0- Using well-known optical trapping tech- 
niques |3|, a ID bosonic experiment would give the two 
beams needed here. 

We consider an initially phase-coherent atomic BEC 
confined to one spatial dimension. (Ref. j2f| gives the 
necessary conditions.) This is then divided into three 
parts. The two outside parts are moved away from the 
central core, and stored for use as local oscillators. The 
central core is coherently converted to a molecular BEC. 
Our theory describes the result of a subsequent dissocia- 
tion of the molecular BEC into two energetic "daughter" 
condensates 0, which interfere with the local oscillators 
to provide the measured quadrature signals. 

The quantum- field-theory effective Hamiltonian de- 
scribing this process is 0, 

H = Ho + H I dz | W*J*< + E 



.x(t) 



(i) 



Here, the atomic and molecular fields are respectively de- 
scribed by the bosonic operators and if? 2, H is the 
kinetic energy, Vi{z) (i = 1,2) are the trapping poten- 
tials (including internal energies), and the Uij are the 
strengths of the one-dimensional intra- and cross-species 
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s-wave interactions. The term x(i) = xo@(ti — t) is re- 
sponsible for coherent conversion of molecules into atom 
pairs, where Xo > and 6{t\ — t) is the Heaviside function 
that turns off the coupling x at t > t\. 

In what follows, we assume that the atom-atom s-wave 
scattering term is negligible. This condition is imposed 
in order to minimize the phase diffusion |2q of the two 
local oscillators required to access the atomic quadra- 
ture correlations. This requires either a low density and 
short interaction times, or else an effective renormalized 
interaction near a magnetic Feshbach resonance which is 
tuned to give a zero effective scattering (Uu = 0). The 
dissociation coupling xo would be caused in the first case 
by a Feshbach sweep, or in the second case [27| by a 
coherent Raman transition with an overall detuning 2A 
Q. This gives an energy mismatch 2HA < between 
the atomic and molecular fields, which is converted into 
kinetic energy [2H\ A| — » 2h 2 k 2 / (2mi)] of atom pairs with 
opposite momenta around fco — ±\/2toi|A|//i, where mi 
is the mass. 

To gain some analytic insight, we will first analyze a 
simple model, beginning with a uniform molecular BEC 
in a coherent state with ID (linear) density rii. The 
condensate extends from —L/2 to L/2, with periodic 
boundary conditions. The dissociation coupling x is 
turned on suddenly, and subsequently assumed to be con- 
stant. Because we are interested in evolution over short 
times, the molecular field depletion is assumed at this 
stage to be negligible so that the amplitude ^2 = \/n2 
(assumed real) can be absorbed into an effective gain 
constant g = xov^2- The dimensionless form of the 
equations is achieved by introducing characteristic time 
and length scales, to — l/ff and e?o = y/^toj (2?ni), and 
transforming to dimensionless time r = t/to, coordinate 
£ = z/do, detuning S — Ato — A/g, and dimensionless 
fields ipi(£,T) = ^i(^d ,Tt )/y/n2- The dimensionless 
initial molecular field density is now scaled to one. 

We expand ipi(Ci T ) m terms of single-mode bosonic 
operators: V>i(£,t) = ^2 q d q (T)e lq ^ /\ft, where q = dgk is 
a dimensionless momentum [/ = L/do, k = {2ir/L)n, n — 
0, ±1 , ±2, ...]. The corresponding Heisenberg equations 
have the solutions a g (r) = A g (r)a g (0) + _B g (r)aL g (0), 
and aigir) = S 9 (r)a 9 (0) + A*(r)aL g (0), where A q (r) = 
cosher) - iX g smh(g q T) /g q , B q (r) = sinh (g q r) /g q , 
X q = q 2 + 5, and g q = (1 — X 2 ) 1 / 2 . The coefficients 
A q and B q satisfy \A q \ 2 — B 2 = 1. The detuning 6 is 
the only parameter that characterizes the dynamics of 
this dimensionless model. For dissociation to proceed, S 
must be negative, which can be achieved by appropriate 
tuning of the frequencies of the Raman lasers. 

In the above solutions, coupling is between opposite 
momentum components only. In quantum optics similar 
solutions have been studied by Reid 0] in the context of 
parametric down-conversion. In that case, the parame- 
ter X q would be identified with an effective phase mis- 



match term, which was set to zero. We note here that, 
unlike photons, the correlated atom pairs are not distin- 
guishable by frequency or polarization, but by different 
momenta or spatial locations. We can now calculate any 
operator moments at time t, given a vacuum initial state 
of the atomic fields. 

We now consider the measurements that must be made 
to demonstrate the EPR paradox. It is readily seen 
that correlations exist between atomic quadratures of the 
beams with opposite momenta. For example, a measure- 
ment of X q — d q + a q allows us to infer, with some error, 

the value of X_ q — a_ g + aL ql and vice versa. The same 

holds for the Y± q = —i(d± q — a± g ) quadratures. This al- 
lows us to define, depending on which beam we measure, 
four inferred variances, 

V mf (X ±q ) = V(X ±q ) [V(X +q ,X_ q )] 2 /V(X Tq ), (2) 

with similar expressions for V m *{Y± q ), where V(a, b) 
n't — (a) (b) . These quadrature correlations can be stud- 
ied using balanced homodyne detection, which mixes the 
signal with a strong local oscillator on the matter-wave 
analog of a 50-50 beam splitter and is a well-known tech- 
nique in quantum optics. Quadratures are measured via 
measuring the density differences, after combining the 
signal and local oscillator [2^ |. 

As an example demonstrating the EPR paradox we 
consider the correlations between the momentum com- 
ponents ±<7o = corresponding to perfect phase 
matching with X q — 0. In this case we obtain 
V inS \X +qa )V inf (Y +qo ) = costT 2 (2T) < I, and the same 
result for V mf (X_ qo )V lnf (Y_ qo ). Since the products of 
the non-inferred variances are bound by the Heisenberg 
uncertainty relation V(X± qo )V(Y± qo ) > 1, this is sim- 
ilar to the EPR paradox. However, as the plane-wave 
momentum components are not localized, this system 
does not allow the required spatial separation of the EPR 
gedanken experiment. 

We now return to the more realistic case described 
by the full Hamiltonian QJ. The dissociated atoms are 
assumed untrapped longitudinally yet confined trans- 
versely, so that they can be treated as a free ID field, 
initially in a vacuum state. The absolute detuning |A| 
must not exceed the trap radial oscillation frequency lu± 
in order to maintain the ID condition [25j with "frozen" 
transverse motion of the dissociated atoms. For com- 
pleteness, the atom-molecule and molecule-molecule scat- 
tering terms are all taken into account. Dimension- 
less interaction couplings are introduced according to 
Ui2 = Ui2y/n2/xoi where = 712(0) is the initial ID 
peak density of the molecular BEC. 

To solve for the resulting quantum dynamics, we use 
stochastic differential equations in the positive-P rep- 
resentation [jfil Ef^ . The essence of the positive-P 
method is in mapping the operator equations of mo- 
tion into stochastic c-number differential equations that 
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can be solved numerically. This requires four inde- 
pendent stochastic fields, ipi and ipf, corresponding to 
the operators ^ and rpl, while i>i(£,t) = u^ip^^ 
and v 2 (£,t) = -tt 22 (l - i 2 Hl) + Y.i u i-2.^t represent 
the effective potentials including the atom-molecule and 
molecule- molecule s-wave interactions. Here, £o = z^/do 
is the dimensionless Thomas-Fermi (TF) radius. We in- 
clude linear losses of atoms and molecules at rates 7$. 
The stochastic variables have a correspondence with nor- 
mally ordered operator moments in the sense of an av- 
erage over a large number of trajectories. The equations 
are 



dt/>i 
dr 

dr 



dipt 



(71 + iS + ivi)ipi + Ktp 2 fpt 
"73), 



\/ nip2 rjx + \J ~iui2ipi^p2/2 (772 

i_ d 2 i>2 
2 



(72 + iv 2 )^2 



\j-iu22^\ Vi + V-^12^1^2/2 (772 



ae 2 



(71 — iS — ivi)ipi + Kipjtp 1 



: 2 <9£ 2 



(72 - iV2)^2 



+ 2 



+ \J iu22^2 2 Vs + \J iuwipti't/ 2 (ve - ivr), (3) 

Here, rjj (j — 1, ...8) are real, independent Gaussian noise 
terms with the correlations ?7j(£, t)%(^', t') — 5jkS(t; — 
£')8(t — t'), and rjj = 0. To numerically integrate these 
equations, we consider that the molecular BEC is ini- 
tially in a coherent state, represented spatially by the 
TF solution. The time duration for the molecule-atom 
conversion is controlled via k(t) = Q(j\ — t), so that 
k(t) — for r > Ti. Once the dissociation stops, we con- 
tinue the evolution of the resulting atomic field in free 
space to allow spatial separation of atoms with positive 
and negative momenta. At the same time we set the 
molecular fields to zero for r > t\ . This models selective 
removal of the molecules by a "blast" pulse with a res- 
onant laser || and is aimed at minimizing the effect of 
atom-molecule scattering and atom losses due to inelastic 
collisions, which can potentially reduce the atom-atom 
correlations. We note that the losses due to inelastic 
collisions are neglected altogether in our model. With 
typical loss rate coefficients of the order of 5 x 10~ 17 
m 3 /s [3] and peak molecular density ~ 10 20 m -3 , their 
disruptive effect should be negligible on submillisecond 
timescales used here. Similar considerations apply to the 
role of losses due to molecule-molecule inelastic collisions 
and three-body losses 0. 

In the nonuniform treatment, due to the mode mixing 
of different momentum components, there is a difference 
in the way the necessary quadratures must be defined 
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Figure 1: (a) Final average atomic density distribution 
ni(£, Tf) (solid line) and the densities of the normalized lo- 
cal oscillator fields |</>+(£)| 2 and |</>_(£)| 2 (dash-dotted and 
dashed lines, with cr+ = 1.45 and cr_ = 1.50, respectively). 
We used a time window of 77 = 70, with dissociation on from 
ro = to ti = 1.65. Other parameter values are S = —25, 
ii22 = 0.2, U12 = 0.1, Co = 238, and j m = 1.2 x 10~ 3 . (b) 
Product of the resulting inferred variances as a function of 
time r for i? = 0.19. The inset shows the same quantity while 
the dissociation is on. 



in comparison to our analytic plane-wave treatment (see 
also Ref. 30]). With the quadratures defined in terms 
of individual Fourier components as above, we obtain no 
inferred violation of the uncertainty principle and hence 
no EPR correlation signature. This is because of the 
assumption of uniform local oscillators, implicitly built 
into this definition of the quadratures. This problem also 
arises in the measurement of pulsed optical squeezing [3lT | 
and is overcome using pulsed local oscillators that are 
mode-matched with the signals. Accordingly, we define 
four mode-matched quadrature operators as 



A±(r) = J de[^(0*i(d-) + ^(0^«,r)], 

Y±(r) = -i [ ^|(0*i(£,T)-<M£)*i(£,T)]. 



(4) 
(5) 



Here, <j>±(£) = |0±(C)l ex P(T*9oC — are t w0 nonuni- 
form local oscillator fields having the same center-of-mass 
momenta, ±qo = ±y/\6\, as the two respective atomic- 
beam signals. These are described classically, with Gaus- 
sian profiles, |</>±(?)| 2 = (27ra 2 = )- 1 / 2 exp(-^ 2 /2cr 2 ; ), with 
/ d£,\<j)±(£,)\ 2 — 1, and are centered at the locations of 
the twin atomic beams at the time of measurement r. 
The phase r) of the local oscillators can be optimized to 
compensate for the molecular mean field phase shift. 

The new quadrature operators have the same com- 
mutation relations as before so that the EPR criterion, 
V mf (X±)V mf (Y±) < 1, and the expressions for the in- 
ferred variances are the same. Experimentally, the fields 
<j)± should be larger than the atomic signals, but they 
have been normalized to one here for convenience. 

Figure QJi shows the final atomic and local oscillator 
densities, obtained from simulation of the full Eqs. 10, 
averaged over 50, 000 stochastic trajectories. In this ex- 
ample |32j| . the average number of atoms in each beam 
after dissociation is 52, with about 10% of them being 
lost during the subsequent free expansion stage. Note 
that, while the two local oscillators should share the 
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same phase, they can have slightly different shapes, or 
atom numbers, without a strong destructive effect on the 
correlations. As we mentioned earlier, the local oscil- 
lators can in principle be prepared by splitting a single 
atomic BEC and then "stored" at spatial locations away 
from the molecular BEC. The zero relative phase can be 
maintained by Feshbach tuning of the magnetic field to 
the zero crossing of the effective atom-atom scattering 
length. In this case, the relative phase drift due to possi- 
ble unequal number of atoms in the two local oscillators 
is minimized. 

Figure shows the product of the inferred variances, 
giving a clear demonstration of the EPR paradox as 
y m f (X + )V mf \Y+) < 1. We note that quadrature cor- 
relations studied here require simultaneous measurement 
of many particles, and that repeated measurements of 
single atomic pairs would not yield the same results. 

To summarize, we have shown that dissociation of a 
molecular BEC into bosonic atoms can provide a simple 
yet robust demonstration of the EPR paradox with mas- 
sive particles. The effects of molecular condensate trap- 
ping and depletion, atom- and molecule-molecule s-wave 
interactions, and possible one-body losses of atoms and 
molecules have all been included in our numerical calcula- 
tions. An experimental realization of our proposal would 
be the first step towards testing fundamental quantum 
mechanics with mesoscopic numbers of massive particles. 

The authors acknowledge support from the Australian 
Research Council, the New Zealand Foundation for Re- 
search, Science and Technology (Grant No. UFRJ0001), 
and thank the authors of the XMDS software 331. 



[1] S. Foiling et ai, Nature 434, 481 (2005). 

[2] M. Greiner et al, Phys. Rev. Lett. 94, 110401 (2005). 

[3] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 

777 (1935); N. Bohr, Phys. Rev. 48, 696 (1935). 
[4] M.D. Reid, Phys. Rev. A 40, 913 (1989); M.D. Reid and 

P.D. Drummond, Phys. Rev. Lett. 60, 2731 (1988). 
[5] Z.Y. Ou et al, Phys. Rev. Lett. 68, 3663 (1992). 
[6] U.V. Poulsen and K. M0lmer, Phys. Rev. A 63, 023604 

(2001). 

[7] K.V. Kheruntsyan and P.D. Drummond, Phys. Rev. A 
66, 031602 (2002); K. V. Kheruntsyan, Phys. Rev. A 71, 
053609 (2005). 

[8] T. Mukaiyama et al, Phys. Rev. Lett. 92, 180402 (2004). 
[9] S. Butt, T. Volz, and G. Rempe, Phys. Rev. A 70, 

03160RR) (2004). 
[10] H. Pu and P. Meystre, Phys. Rev. Lett. 85, 3987 (2000); 

L.-M. Duan et ai, ibid. 3991 (2000). 
[11] R.H. Wynar et al, Science 287, 1016 (2000). 
[12] E.A. Donley et ai, Nature 417, 529 (2002); J. Herbig et 

ai, Science 301, 1510 (2003); K. Xu. et al, Phys. Rev. 

Lett. 91, 210402 (2003); S. Diirr et ai, ibid. 92, 020406 

(2004). 

[13] P. D. Drummond, K. V. Kheruntsyan, and H. He, Phys. 
Rev. Lett. 81, 3055 (1998); K.V. Kheruntsyan and P.D. 



Drummond, Phys. Rev. A 61, 063816 (2000). 
[14] D.J. Heinzen, et al. Phys. Rev. Lett. 84, 5029 (2000); 

P.D. Drummond, K.V. Kheruntsyan, D.J. Heinzen and 

R.H. Wynar, Phys. Rev. A 65, 063619 (2002). 
[15] E. Timmermans et al, Phys. Rev. Lett. 83, 2691 (1999); 

E. Timmermans et ai, Phys. Rep. 315, 199 (1999). 
[16] J. Javanainen and M. Mackie, Phys. Rev. A 59, R3186 

(1999); M. Mackie, R. Kowalski and J. Javanainen, Phys. 

Rev. Lett. 84, 3803 (2000). 
[17] V. A. Yurovsky et al, Phys. Rev. A 60, R765 (1999). 
[18] F.A. van Abeelen and B.J. Verhaar, Phys. Rev. Lett. 83, 

1550 (1999). 

[19] M. Holland, J. Park, and R. Walser, Phys. Rev. Lett. 86, 

1915 (2001); S.J.J.M.F. Kokkelmans and M.J. Holland, 

ibid. 89, 180401 (2002). 
[20] J.J. Hope and M.K. Olsen, Phys. Rev. Lett. 86, 3220 

(2001); M.K. Olsen and L.I. Plimak, Phys. Rev. A 68, 

031603 (2003); 

[21] K. Goral, M. Gajda, and K. Rzazewski, Phys. Rev. Lett. 

86, 1397 (2001). 
[22] T. Kohler, T. Gasenzer, and K. Burnett, Phys. Rev. A 

67, 013601 (2003). 
[23] R.A. Duine and H.T.C. Stoof, Phys. Rev. Lett. 91, 

150405 (2003); Phys. Rep. 396, 115 (2004). 
[24] L. Tolra et al., Phys. Rev. Lett. 92, 190401 (2004); T. 

Kinoshita et al., Science 305, 1125 (2004). 
[25] K.V. Kheruntsyan et al., Phys. Rev. Lett. 91, 040403 

(2003); Phys. Rev. A 71, 053615 (2005). 
[26] M.J. Steel, M.K. Olsen, et ai, Phys. Rev. A 58, 4824 

(1998). 

[27] In a different context, the strategy of combining a Fesh- 
bach resonance and Raman transitions for atom-molecule 
conversion has been considered by S.J.J.M.F. Kokkel- 
mans et al., Phys. Rev. A 63, 031601 (2001) and M. 
Mackie, ibid. 66, 043613 (2002). 

[28] H.-A. Bachor, A guide to experiments in quantum optics, 
(Wiley- VCH, Weinheim, 1998). 

[29] P.D. Drummond and CW. Gardiner, J. Phys. A 13, 2353 
(1980); P.D. Drummond and JF. Corney, Phys. Rev. A 
60, 2661(R) (1999). 

[30] V.A. Yurovsky and A. Ben-Reuven, Phys. Rev. A 67, 
043611 (2003). 

[31] M.G. Raymer, P.D. Drummond, and S.J. Carter, Opt. 
Lett. 16, 1189 (1991). 

[32] Considering a 87 Rb2 experiment as a guideline, these pa- 
rameters relate to a molecular BEC confined by a har- 
monic trap with an axial (radial) frequency uj z /2tv = 1 
Hz (cu±/2n — 12 kHz). We chose the molecule-molecule 
scattering length 022 = 1 nm, giving the ID coupling 
U22 — 2ha22/m2l 2 L — 1.5 x 10~ 4 m s -1 , where 7712 = 2mi 
and lx = \J h/m2LJj_. In addition, we take the atom- 
molecule s-wave interaction strength U12 = 0.75 x 10~ 4 
m s - ; the relatively small values of these couplings 
are chosen here for making the computation feasible. 
The initial number of molecules is N2 = 380, at ID 
peak density of n 2 (0) = nl ± ri£ D) (0) = 3.1 x 10 6 m" 1 
{nf D) (0) = 2 x 10 20 m" 3 ), while the TF radius is 
zo — \J 2hU22ri2 (0) / {miU3^) = 93 fim . We take the de- 
tuning A = —6.1 x 10 4 s~\ and the atom-molecule cou- 
pling x (3D) = 3 .4 x 1 0~ 7 m 3/2 s" 1 , so that the ID value 
is % — X (3 ° ) /V / 4^± = 1-4 m 1/2 s" 1 . The characteristic 
time scale is to = 0.41 ms, so that the duration of dis- 
sociation (n = 1.65) corresponds to ii = 0.68 ms, while 



5 

the total time window (77 = 70) is tf = 29 ms. Hope, www.xmds.org. 

[33] P.T. Cochrane, G. Collecutt, P.D. Drummond, and J.J. 



